An and Lin BMC Plant Biology 201 1, 11:105 
http://www.biomedcentral.eom/1 471-2229/1 1 /1 05 



Plant Biology 



RESEARCH ARTICLE Open Access 



Transcriptional regulatory programs underlying 
barley germination and regulatory functions of 
Gibberellin and abscisic acid 

Yong-Qiang An 1 " and Li Lin 2 



Abstract 

Background: Seed germination is a complex multi-stage developmental process, and mainly accomplished 
through concerted activities of many gene products and biological pathways that are often subjected to strict 
developmental regulation. Gibberellins (GA) and abscisic acid (ABA) are two key phytohormones regulating seed 
germination and seedling growth. However, transcriptional regulatory networks underlying seed germination and 
its associated biological pathways are largely unknown. 

Results: The studies examined transcriptomes of barley representing six distinct and well characterized 
germination stages and revealed that the transcriptional regulatory program underlying barley germination was 
composed of early, late, and post-germination phases. Each phase was accompanied with transcriptional up- 
regulation of distinct biological pathways. Cell wall synthesis and regulatory components including transcription 
factors, signaling and post-translational modification components were specifically and transiently up-regulated in 
early germination phase while histone families and many metabolic pathways were up-regulated in late 
germination phase. Photosynthesis and seed reserve mobilization pathways were up-regulated in post-germination 
phase. However, stress related pathways and seed storage proteins were suppressed through the entire course of 
germination. A set of genes were transiently up-regulated within three hours of imbibition, and might play roles in 
initiating biological pathways involved in seed germination. However, highly abundant transcripts in dry barley and 
Arobidopsis seeds were significantly conserved. Comparison with transcriptomes of barley aleurone in response to 
GA and ABA identified three sets of germination responsive genes that were regulated coordinately by GA, 
antagonistically by ABA, and coordinately by GA but antagonistically by ABA. Major CHO metabolism, cell wall 
degradation and protein degradation pathways were up-regulated by both GA and seed germination. Those genes 
and metabolic pathways are likely to be important parts of transcriptional regulatory networks underlying GA and 
ABA regulation of seed germination and seedling growth. 

Conclusions: The studies developed a model depicting transcriptional regulatory programs underlying barley 
germination and GA and ABA regulation of germination at gene, pathway and systems levels, and established a 
standard transcriptome reference for further integration with various -omics and biological data to illustrate 
biological networks underlying seed germination. The studies also generated a great amount of systems biological 
evidence for previously proposed hypotheses, and developed a number of new hypotheses on transcriptional 
regulation of seed germination for further experimental validation. 
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Background 

Seed germination is a complex multi-stage developmen- 
tal process important to plant development, plant evolu- 
tion, and agricultural production. Strictly defined, 
germination begins with the uptake of water by dry 
quiescent seeds and ends with the visible emergence of 
an embryo tissue from its surrounding tissues. However, 
in many scientific literatures and agronomic research, 
seed germination often broadly includes early seedling 
growth, a process which ends with the start of auto- 
trophic growth or the emergence of seedling from soil 
[1]. Seed germination is accompanied with many distinct 
metabolic, cellular and physiological changes. For exam- 
ple, upon imbibition, the dry quiescent seeds take up 
water and rapidly resume many fundamental metabolic 
activities such as respiration, RNA and protein synthesis 
machinery, as well many enzyme activities using surviv- 
ing structures and components in the desiccated cells. 
Meanwhile, dry seeds gradually lose stress tolerances, 
such as desiccation tolerance, over the course of seed 
germination. These combined biological activities trans- 
form a dehydrated and resting embryo with an almost 
undetectable metabolism into one with vigorous meta- 
bolism calumniating in growth [2,3]. 

GA and ABA are two key phytohormones regulating 
seed germination and seedling growth. It is believed that 
GA and ABA play antagonistic roles in regulating seed 
germination and their ratios govern the maturation ver- 
sus germination pathways that embryos will take after 
they complete rudimentary organogenesis [4,5]. It was 
proposed that GA enhances seed germination and seed- 
ling growth. Maturing maize embryos require GA for 
germination in culture. Treating maize embryos with GA 
synthesis inhibitors also decrease both the rate of germi- 
nation and the fraction of embryos that germinate [4]. 
Treatments that promote Arabidopsis germination, such 
as cold and light, are often correlated with an increase in 
endogenous GA [6]. It has been showed that GA-defi- 
cient Arabidopsis and tomato mutants are impaired in 
seed germination [7,8]. It is proposed that a conserved 
DELLA protein negatively mediates GA regulation of 
seed germination and seedling growth [9-13]. However, 
the biological networks underlying GA regulation of seed 
germination and seedling growth are largely unknown. In 
germinating cereal grains, GA is primarily synthesized in 
the embryo and is then relocated to aleurone tissues 
where it induces synthesis of hydrolytic enzymes. The 
hydrolytic enzymes are further secreted into starchy 
endosperm to mobilize seed storage reserve to provide 
nutrients and energy for embryo growth and differentia- 
tion before an autotrophic phase is fully established. It is 
believed that GA induction of hydrolytic activities mainly 
occurs in the post-germination phase to support seedling 
growth [14]. However, the requirement of GA in early 



barley germination remains to be determined. In con- 
trast, ABA content increases dramatically in most plant 
species during seed maturation, and induces the produc- 
tion of seed storage and desiccation tolerant proteins to 
prepare the seeds for undergoing desiccation and to pro- 
duce energy and nutrient reserve for later seed germina- 
tion [15-18]. ABA also suppresses expression of many 
hydrolytic enzyme genes to prevent viviparous germina- 
tion [19,20]. Recent evidence suggests that other phyto- 
hormones, including auxin and ethylene, play roles in 
regulating seed germination [21,22]. 

Barley germination and seedling growth have been 
investigated extensively due to its importance in barley 
agriculture and the brewing industry [23]. A wealthy 
amount of diverse biological data from barley germina- 
tion has been accumulated [24]. In addition, barley 
aleurone from germinating barley grains has been estab- 
lished as a model system to study the mode of action on 
GA and ABA response pathways and their regulatory 
functions in barley seed germination [25]. Recently, Bar- 
ley Genome GeneChips containing approximately 
22,700 genes were used in examining the transcriptome 
of barley aleurone in response to GA, ABA, and the 
inactivation of SLN1 proteins. The analysis identified 
1328 GA and 206 ABA responsive genes and revealed 
that transcriptomes of barley aleurone respond antago- 
nistically to GA and ABA treatments [20]. Loss-of- func- 
tion of the DELLA protein, SLN1, activates barley 
aleurone transcriptomic programs in response to GA 
[26]. A great number of studies examined transcrip- 
tomes of Arabidopsis germinating seeds and tissues to 
study seed germination in response to developmental 
regulation, genetic variation, and environmental signals 
at a system level [27-30] . Transcriptomes of various ger- 
minating tissues in barley have been determined using a 
variety of transcript profiling technologies [31-35]. Bar- 
rero et.al compared the transcriptomes of coleorhiza 
and roots from dormant and after-ripened barley 
embryos at 8 and 18 hours after imbibition and charac- 
terized the dormancy related transcriptomic changes. 
Screenvasulu et. al performed a transcriptome analysis 
of endosperm and embryo at barely grain maturation, 
desiccation, and early seedling growth stages, and 
revealed a smooth transition in the transcriptional pro- 
gram between late seed maturation and early seedling 
growth within embryo tissues [33]. However, the 
research mainly focuses on post-germination processes. 
No germinating barley tissues prior to emergence of 
coleorhiza from grains were examined. It has been well 
demonstrated that the activities of many germination 
related gene products and biological pathways are sub- 
ject to strict developmental regulation. However, an in- 
depth and comprehensive transcriptomic characteriza- 
tion of germinating barley representing distinct and well 



An and Lin BMC Plant Biology 201 1, 11:105 
http://www.biomedcentral.eom/1 471-2229/1 1 /1 05 



Page 3 of 24 



defined developmental stages over the entire course of 
seed germination are not available. To fill the gap, the 
studies carefully examined several physiological and 
morphological characteristics of barley germination and 
seedling growth, and selected six distinct developmental 
stages that represent the entire process of barley germi- 
nation from grain imbibition to early seedling growth 
for transcriptome analysis. Extensive bioinformatic ana- 
lysis of the dynamic transcriptomic data delineates the 
transcriptional regulatory program underlying barley 
germination and seedling growth at gene, pathway and 
systems levels. 

Results and Discussion 

Distinct Physiological and Developmental Stages of 
Barley Germination 

One of the experimental objectives is to determine 
dynamic changes in transcriptomes of barley over the 
course of seed germination, and to further illustrate the 
transcriptional regulatory program underlying barley 
germination and its associated biological pathways. 
However, expression of germination important genes 
and biological pathways is often subjected to strict 
developmental regulation over the course of seed germi- 
nation. It is crucial to examine transcriptome of barley 
representing well defined and distinct physiological 
stages. Morphology, water up-take, amylase activity and 
loss of seed desiccation resistance are important charac- 
teristics of germinating barley [2]. Having examined 
these characteristics of germinating barley over the 
course of seed germination, six developmental stages 
representing distinct physiology of barley over the 
course of seed germination were selected and referred 
as SO to S5 stages. 

Figure 1A shows the morphology of germinating bar- 
ley at each developmental stage and time typically taken 
for dry mature grains to reach the given stage. For 
example, S3 stage marked the end point of germination 
process, and can be easily identified by the visible 
coleorhiza emergence from the grains, which typically 
occurred at 18 hours of germination. No morphological 
changes were observed for the germinating grains prior 
to the S3 stage. Germinating grains at time points of 
1/6 and 1/2 of the time typically taken for coleorhiza to 
emerge from germinating grains were referred to SI and 
S2 stages to represent the early stages of germination. 
Figure IB shows that water uptake of germinating barley 
had three phases over the course of seed germination as 
previously described [1]. The water content of germinat- 
ing barley rapidly increased between SO and SI stages at 
a rate of 7.1% per hour. However, the water uptake slo- 
wed down dramatically after the SI stage to the lowest 
rate of 1.2% per hour between S2 and S3 stages. Follow- 
ing the low-point, water uptake gradually increased to a 



higher rate of 5.2% per hour between stages S4 and S5 
(Figure IB). Alpha-amylase activity exemplifies the 
mobilization of starch storage reserves over the course 
of seed germination. While there was little change in 
alpha-amylase activity until the S3 stage, a dramatic 
activity increase occurred between stages S3 and S4, fol- 
lowing the emergence of coleorhiza from germinating 
grains (Figure 1C). Mature dry seeds are highly resistant 
to desiccation and many other abiotic and biotic stres- 
ses. Over the process of grain germination, grains gradu- 
ally lose desiccation resistance and those stress 
tolerances. No significant change in desiccation resis- 
tance was observed for the germinating grains until the 
S2 stage (Figure ID). However, the survival rate of dehy- 
drated grains at S3 stage dropped to 14.7% (Figure ID). 
No germinating barley at S4 could be revived after dehy- 
dration. Barley grains completely lost their desiccation 
resistance over the period after grains finished their ger- 
mination at S3 stage and before the production of amy- 
lase increases significantly at S4 stage. Thus, the six well 
characterized developmental stages defined above should 
cover the entire spectrum of physiological changes in 
barley from initial grain imbibition to early seedling 
growth. The germinating barley at each of the six stages 
should provide a representative and distinct physiologi- 
cal and developmental stage, and can be accurately and 
easily identified based on the relative timing of germina- 
tion and morphology of seedlings. 

It is a great challenge to accurately define and identify 
physiology and developmental stages of germinating 
seeds. Time points post imbibition are widely used to 
define the developmental stages of seed germination. 
However, seed germination rates are significantly 
affected by genotypes, physiology of dry seeds and ger- 
mination environments. Individual dry seeds from the 
same harvest do not always germinate uniformly due to 
heterogeneity of seed maturity [2]. Although a large 
amount of diverse biological data and results have been 
reported in the previous germination studies, it faces a 
great difficulty to compare or integrate those data 
because the developmental stages of the germinating tis- 
sues used in most of those studies are not well defined. 
The relative timing of germination and the morphology 
of seedling described above could be used as an accurate 
and facile approach to identify germinating grains and 
seedlings equivalent to each of the six developmental 
stages and control variation of germination rates caused 
by those factors in other cereal species. 

A Transcriptomic Switch Correlated to the Morphological 
and Physiological Transition from Seed Germination to 
Seedling Growth 

Affymetrix Barley Genome GeneChip Arrays containing 
22,792 probe sets [36] were used to examine the 
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Figure 1 Morphology and Physiology of Germinating Barley at Each Developmental Stage. Figure 1A. Morphology and time points of 
germinating barley at each developmental stageThe morphology of germinating grains and seedlings at each developmental stage and the 
typical time taken for dry mature grains to reach each stage are shown. The relative times in reference to the time taken to reach S3 stage are 
indicated in the parenthesis. S3: Coleorhiza emerging from grains at 18 hours of germination. S4: the rootlet length is half that of its grain at 33 
hours of germination. S5: the shoot is 3 times as long as its grain at 71 hours of germination. Figure 1B. Water content of germinating barley at 
each developmental stageThe fresh and dried weights of 10 germinating grains or seedlings at each developmental stages defined as in the 
Figure 1 A were measured. The water content in germinating barley at each stage is indicated as Y axis as a percentage of the dry weight. The 
representative time point of germination at each stage is indicated as X axis. Standard derivations of the three replications are indicated as error 
bars. Stages are marked. Figure 1C Alpha-amylase activity in germinating barley at each developmental stageThe X axis indicates the 
development stages. The average amount of maltose in umole produced per gram of fresh examined tissues (U/gfw) and the standard 
derivation of three replications are indicated on Y axis. Figure 1D. Desiccation resistance of germinating barley at each developmental stageThe 
germinating barley at each developmental stage were dehydrated, and then re-germinated. The percentage of the dehydrated germinating 
barley that could revive to their growth was defined as survival rate and indicated as Y axis to measure desiccation resistance of the germinating 
barley. The X axis indicates developmental stages. Standard derivations of the three replications are indicated as error bars. 



transcriptomes of germinating barley at the SO, SI, S2, 
S3, S4 and S5 stages. Three GeneChip assay replications 
each from an independent germination experiment were 
conducted for each developmental stage to control bio- 
logical and technical variation. Figure 2 summarizes the 
number of mRNA species accumulating at a detectable 
level at each developmental stage. Over 50% of the 
22,792 examined transcript species were accumulated at 



detectable levels in the dry grain. However, mRNA com- 
plexity increased over the course of germination except 
for a decrease from the SI to S2 stages. The most dra- 
matic increase in mRNA complexity occurred between 
the S2 and S3 stages. It is consistent with the previous 
reports that approximately 50% of examined transcripts 
are accumulated at a detectable level and encode all 
functional categories of proteins in dry seeds of 
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SO SI S2 S3 S4 S5 

Figure 2 Total Number of Detectable Transcript Species in Germinating Barley At Each Developmental Stage A probe set with present 
calls at P value cutoff of 0.05 in two of three replicates at each developmental stage is assigned as a transcript expressed at a detectable level. 
The number of probe-sets with "present" calls at each developmental stage is indicated as the Y-axis. The X-axis indicates the developmental 
stages. 



divergent plant species [27,28,33]. Although a number of 
transcripts encoding seed maturation specific proteins 
such as seed storage proteins have been demonstrated 
to degrade over the course of germination [31], it is 
likely that many of them are preserved and continue to 
function through seed germination, at least through the 
early seed germination. It was shown that germination 
of Arabidopsis seeds can be blocked by a translational 
inhibitor, cycloheximide, but not by a RNA polymerase 
II inhibitor, alpha-amanitin [37,38]. Thus, it is likely that 
the potential of germination is largely programmed in 
the seed developmental process. 

GC-RMA algorithm was used to convert probe level data 
to expression measurement in the microarray experiments 
[39]. One-way ANOVA analysis identified 6157 genes 
whose transcript accumulation changed significantly over 
the process of barley germination with a False Discovery 
Rate (FDR) of 0.05 (See Additional file 1). 5382 genes were 
differentially regulated between SO stage and any other 
developmental stage. Of the 5382 genes, 4493 genes (84%) 
showed more than a three-fold change, indicating that 
most of the differentially regulated genes changed dramati- 
cally over the course of seed germination (See Figure 3 and 
See Additional file 2). Of the 4493 genes, 2,816 genes were 
up-regulated while 1,688 genes were down-regulated. 
There were 63% more up-regulated genes than down-regu- 
lated genes. This is consistent with increasing complexity 
of mRNA species and biological processes in germinating 
grains over the course of germination. 
Comparison of any two adjacent developmental stages 
shows that 3267 genes changed significantly in their 



mRNA accumulation (Figure 3 and Additional file 2). 
2390 of the 3267 genes had more than a three-fold 
mRNA accumulation changes. Interestingly, 1295 genes 
showed differential accumulation of their mRNAs within 
the 9 hours of germination between the S2 and S3 
stages while only 310 genes were differentially expressed 
within following 15 hours between S3 and S4 stages. 
Both the biggest increase in mRNA complexity and the 
highest number of differentially regulated genes were 
observed between the S2 and S3 stages. Thus, a dra- 
matic transcriptional program switch occurred between 
the two developmental stages and co-occurred with 
morphological emergence of coleorhiza, dramatic 
decrease in desiccation resistance, initiation of enzy- 
matic alpha amylase activity increase, and the slowest 
water uptake over the course of germination. The tran- 
scriptional reprogram switch is likely to play a key regu- 
latory role in transforming barley grains from 
germination to seedling growth. 

The majority of differentially regulated genes between 
adjacent developmental stages showed more than three- 
fold changes in their transcript accumulations. To focus 
on the genes that are more likely to have functional sig- 
nificance in seed germination, the following analysis and 
description are only limited to the genes with more than 
three-fold changes unless specified otherwise. 

Conservation and Divergence of Highly Abundant 
Transcripts in Barley and Arabidopsis Dry Seeds 

Table 1 lists top 100 barley probe-sets that have the 
highest signal intensity in barley dry grains. The probe- 
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Figure 3 Summary of Differentially Regulated Genes over the Course of Germination. Differentially regulated genes were identified by 
either comparison of adjacent stages (left 6 bars) or comparing each stage to SO stage (right 6 bars). The numbers of up- or down-regulated 
genes with statistical significance but less than three-fold change in their mRNA accumulation are indicated in pink or light blue. The numbers 
of up- or down-regulated genes with both statistical significance and more than a three-fold change in their mRNA accumulation are indicated 
in red or blue. The total number of genes differentially regulated between any adjacent stages and between any stage and SO stage are marked 
as "Any adjacent" and "any stage/SO". 



sets accounted for less than 1% of transcripts detectable 
in dry grains, and should represent highly abundant 
transcripts stored in the dry seeds. Those transcripts 
mainly encode proteins related to nutrient reservoir, 
stress tolerance, protein biosynthesis, glycolysis, lipid 
metabolism, oxidoreduction, and metal binding. A num- 
ber of transcripts encoding proteins with unknown func- 
tions or related to other diverse functionalities were also 
found in the top 100 barley probe-sets. 

Transcriptomes of Arabidopsis dry seeds have also 
been extensively characterized [28,29,38]. Nakabayashi 
et. al identified 484 highly abundant transcripts in non- 
dormant Arabidopsis dry seeds. Comparing the two sets 
of Arabidopsis and barley genes identified 35 pairs of 
putative barley-Arabidopsis orthologs that are highly 
abundant in both barley and Arabidopsis dry seeds. In 
addition, ten pairs of homologous bzucley-Arabidopsis 
genes with e-value less than -10 have been identified in 
those functional categories. Those transcripts were 
found in all functional groups of highly abundant barley 
transcripts except for the groups of unknown functions 
and glycolysis pathways. Many of their encoded proteins 
and pathways have been previously reported to be highly 
accumulated in the dry seeds, and suggested to be 
involved in seed maturation and germination [40]. For 
example, increasing evidence indicates that germination 
of seeds is accompanied by extensive changes in the 
redox state of proteins [41,42]. Translation of dry seed 



stored transcripts is required for seed germination 
[37,38]. Monocot-dicot divergence occurred approxi- 
mately 200 million years ago [43]. Gene expression pat- 
terns change quickly if they have no functional 
constrains [44-47]. Preserving high accumulation of 
those ancient gene transcripts and pathway transcripts 
in both barley and Arabidopsis dry seeds from their 
ancestor after 200 million years of independent evolu- 
tion strongly suggests that those transcripts and path- 
ways are functionally important to germination, and 
may contribute to the biological characteristics of germi- 
nation shared by barley and Arabidopsis, Although bar- 
ley and Arabidopsis have evolved as two distinct types 
of starchy and oil seed plants respectively over the 200 
million years, it is likely that transcriptional programs 
and molecular mechanism underlying seed germination 
are highly conserved, particularly in biological pathways 
such as stress tolerance, nutrient reservoir and protein 
translation. Interestingly, two pairs of oleosin ortholo- 
gous transcripts are highly accumulated in not only Ara- 
bidopsis seeds but also barley seeds. Oleosin is a highly 
accumulated protein in oil bodies that mainly stores 
triacylglycerol (TAG) as major reserve in mature seeds 
to provide energy for seed germination and seedling 
growth [48]. It is believed that oleosin plays important 
regulatory roles in oil body stabilization and size [49,50]. 
It was observed that oleosin proteins are highly abun- 
dant in oil bodies from Arabidopsis and Brassica seeds 
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Table 1 Comparison of Highly Abundant Seed Stored Transcripts in Barley and Arabidopsis (Continued) 
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1: rank of barley transcript abundance in dry seeds. 
2: rank of Arabidopsis transcript abundance in dry seeds 

3. Presence of Arabidopsis genes orthologous to the barley genes: Yes (Y) and No (N) 
*: Arabidopsis genes homologous to the barley genes at e-value less than -10. 
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[51,52]. Oil bodies and expression of oleosin have also 
been observed in barley embryo and aleurone tissues 
[53]. Although barley and Arabidopsis evolve to use 
starch and oil as major storage reserve respectively to 
support seed germination and seedling growth, it seems 
that barley still preserve high accumulation of oleosin in 
seeds. It will be interesting in understanding their biolo- 
gical funcitons. 

A significant number of the highly abundant barley 
seed transcripts have no orthologous genes in Arabidop- 
sis; or their orthologs or strong homologs do not highly 
accumulate in Arabidopsis dry seeds. Some of the barley 
transcripts encode hordein proteins in nutrient reser- 
voir, glycolysis pathway enzymes, proteins with 
unknown functions, and a number of proteins with 
other functions. Interestingly, fructose-bisphosphate 
aldolase and enolase transcripts in the glycolysis path- 
ways are highly accumulated in the barley grains, but 
none of their Arabidopsis orthologs and strong homo- 
logs highly accumulates in Arabidopsis dry seeds. Thus, 
specific high accumulation of the glycolysis enzyme 
transcripts in starch barley dry grains suggest that barley 
has evolved an unique regulatory pathway to quickly 
activate glycolysis upon imbibition to support early 
energy-demanding biological process. It raises possibili- 
ties that those barley genes and/or their high accumula- 
tion patterns in dry grains have diverged from their 
Arabidopsis orthologs after monocot-dicot occurred, 
and contribute to characteristics of barley seeds distinct 
from that of Arabidopsis. The comparative studies on 
the highly abundant barley and Arabidopsis transcripts 
should provide insight into molecular mechanism 
underlying conserved and divergent characteristics of 
barley and Arabidopsis germination. 

The Early and Transient Regulation of Barley Germination 

Transcriptional changes occurred as early as in the first 
three hours of germination. Forty -seven genes were dif- 
ferentially regulated between SO and SI stages. Twenty- 
five of these genes had more than 3 fold increases in 
their mRNA accumulation. Ten of the 25 up-regulated 
genes reached the highest expression level at SI stage 
(Figure 4A), and then gradually dropped to the levels of 
mature grains at S3 stage. This group of genes encoded 
two zinc finger proteins, one Avr9/Cf-9 rapidly elicited 
protein, one DRE-binding protein, one arabinogalactan- 
like protein, two glutaredoxin and three proteins with 
unknown functions. The accumulation of the other 15 
gene transcripts increased at SI stage and reached the 
maximum level at S2 stage (Figure 4B). Those genes 
encoded WRKY family transcription factors, DnaJ-like 
proteins, an Avr9/Cf-9 rapidly elicited protein, a beta- 
glucan elicitor receptor, AAA-type ATPase, serine/ 
threonine phosphatase 2C, ARM repeat protein, 



oxysterol-binding protein-like protein and proteins with 
unknown functions. Interestingly, the transcript accu- 
mulation for the majority of these genes also dropped to 
the levels of mature seeds at S3 stage. Many of the early 
induced genes encoded transcription factors and recep- 
tor proteins. Early differential expression of genes in 
response to GA has been successfully used as a criterion 
to identify regulatory genes in the GA response pathway 
in Arabidopsis [54]. It was shown that transcriptional 
changes can be detected in 15 minute of imbibition in 
Arabidopsis. However, much lower number of genes are 
up-regulated than down- regulated in Arabidopsis seeds 
within the first hour of imbibition. Only four transcripts 
are up-regulated while eighty-three transcripts are 
down-regulated within the first hour of imbibition [29]. 
Such a transient and early induction accumulation pat- 
tern was also observed in the rice germination. A cluster 
of rice transcripts are up-regulated at early stage of ger- 
mination, and reach its peak in their abundance at 1 or 
3 hours after imbibition, and then decreased to low 
levels again at 12 hours after imbibition [55]. The early 
and transient induction of those genes during seed ger- 
mination raises a possibility that the genes could poten- 
tially play regulatory roles in initiating transcriptional 
regulatory cascades and signaling transduction pathways 
underlying barley germination. Some of the genes 
encode transcription factors and regulatory components 
in signaling pathways that are potentially related to seed 
germination. For examples, the probe sets, Contig 
9265_at, encodes a serine/threonine phosphatase type 
2C (PP2C), and was up-regulated by 3.5 and 4 folds at 
SI and S2 stages respectively. It has been shown that 
serine/threonine protein phosphatase 2Cs (PP2C) can 
suppress ABA signaling pathways in Arabidopsis [56]. 
The loss-of- function of Arabidopsis ABI1 and ABI2 that 
encode protein phosphatase 2Cs increases seed dor- 
mancy and enhances responsiveness to ABA [57]. In 
addition, it was reported that two rice WRKY genes 
could repress ABA induction of the HVA22 promoter 
[58]. It is well documented that ABA promotes seed 
dormancy and inhibits seed germination and seedling 
growth [59,60]. The early and transient transcriptional 
up-regulation of negative regulators in ABA signaling 
pathways suggests that the induced accumulation of the 
mRNA species might suppress ABA function at the 
early stage of seed germination to promote seed 
germination. 

Three Distinct Phases of Transcriptional Regulatory 
Program Underlying Barley Germination 

Hierarchical clustering of all examined stages based on 
the normalized mRNA accumulation of the 6157 differ- 
entially regulated genes revealed that the six develop- 
mental stages were further clustered into two groups 
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Figure 4 Expression Patterns of the Genes Up-Regulated in The first Three hours of Germination. The expression pattern of each gene 
with peak expression at S1 stage is shown in 4A and those with peak expression at S2 are shown in 4B. The signal intensity of each gene is 
expressed on the Y-axis. The probe-set ID for each gene is shown. 
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with the threshold distance of 174 (Figure 5). The devel- 
opmental stages of SO, SI, and S2 were clustered into 
one group while S3, S4 and S5 were grouped into 
another. Although S2 and S3 stages were development 
tally adjacent to each other with only 9 hours of inter- 
val, they were clustered in two separate groups, which is 
consistent to the dramatic increase in mRNA complexity 
and the highest number of differentially regulated genes 
observed between the two adjacent developmental 
stages. The clustering data further supports that a dra- 
matic transcriptional program switch occurred between 
the two developmental stages. It also revealed that the 
transcriptional regulatory program underlying germina- 
tion was composed of three distinct transcriptional 
phases of germination. The three distinct phases of the 
transcriptional regulatory program are well correlated to 
the three phases of water up-takes of seed germination, 
and are referred to early (from SO to S2), late (S2 to S3) 
and post- (S3 to S5) germination phases. 

A total of 730, 1295, and 1394 genes changed signifi- 
cantly for more than three folds in their transcript accu- 
mulation during early, late and post-germination phases 
respectively (Additional file 3). Those genes were named 
as early germination, late germination and post-germi- 
nation regulated genes. Figure 6 shows the Venn dia- 
gram of the early, late and post-germination genes with 
426 genes specifically and differentially regulated during 
the early germination phase, 792 genes during the late 
germination phase, and 1051 genes during the post-ger- 
mination phase (See Additional file 3). It is likely that 
those genes are responsible for the biological changes 
specifically occurring in their corresponding germination 



phase. 42 genes were differentially regulated by all three 
developmental phases (Additional file 3), while 51 genes 
by both early and post-germination phases, 250 by late 
and post-germination phases and 211 genes by early and 
late germination phases. 

Biological Pathways Differentially Regulated in Each 
Germination Phase 

Identification of a large number of differentially regu- 
lated genes in each germination phase enables the stu- 
dies to identify biological pathways and functional 
groups that are transcriptionally up-regulated or down- 
regulated at a systems level using MapMan and Page- 
Man tools. In the analysis, the probe-sets on Barley 
Genome GeneChip were assigned into 35 major func- 
tional bins based on their molecular functionalities such 
as metabolic pathways, signaling pathways and gene 
families. Each bin is further divided into sub-functional 
bins [33]. PageMan and MapMan software tools were 
used to determine the statistic probability of over- or 
under-representation of the early, late or post- germina- 
tion regulated genes in each bin and sub-bin [61,62]. 
The representation analysis revealed that 18 bins and 
138 sub-functional bins were up-regulated or down- 
regulated by at the least one of the three germination 
phases with over-representation Z-value of more than 1 
(Figure 7 and 8). Having only considered the pathways 
with statistical significance, the systems and over-repre- 
sentation analysis offers a more effective approach to 
discover biological pathways/processes that are tran- 
scriptionally regulated during seed germination and 
potentially important to seed germination. 
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SO SI 
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Late germination 




Figure 5 Three Phases of the Transcriptional Regulatory Program Underlying Barley Germination. The normalized mRNA intensities of 
the 6157 differentially regulated genes at the six developmental stages were subject to hierarchical cluster analysis. The six developmental 
stages were clustered into two groups with a threshold distance of 174. The transcriptional regulatory program underlying barley germination 
process was divided into three germination phases, early, late and post-germination phases. 
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Figure 6 Venn Diagram of the Early, Late, and Post-germination Regulated Genes. The genes differentially regulated by early (blue), late 
(green), and post- (orange) germination phases were compared and displayed as a Venn diagram. The numbers of differentially regulated genes 
specific to one, two or all processes are indicated in the Venn diagram. 



Up-regulation of Regulatory Components and Cell Wall 
Metabolism in Early Germination Phase 

The representation analysis of the 730 early germination 
regulated genes revealed that bin and sub-bins with reg- 
ulatory and signaling functionalities were preferentially 
up-regulated in the early germination phase (Figure 7 
and 8). The early germination up-regulated genes were 
over-represented in the bins of Cell Wall, RNA and Sig- 
naling. Although the RNA bin is composed of three 
sub-functional bins, RNA Processing, RNA Transcrip- 
tion, and Regulation of Transcription [61], only the Reg- 
ulation of Transcription sub-bin was preferentially up- 
regulated in the early germination phase. All up-regu- 
lated sub-bins in Regulation of Transcription were tran- 
scription factor families including Helix-Loop-Helix, 
C2C2(Zn)GATA, GRAS, NAC domain, and WRKY 
transcription factor families. In the Signaling bin, the 
sub-functional bins of Signaling in Sugar and Nutrient 
Physiology and Receptor Kinases were up-regulated. 
Sub-functional bin of Protein Post-Translational Modifi- 
cation, which often plays regulatory functions in activ- 
ities of proteins, was also preferentially up-regulated in 
early germination phase. The regulatory and signaling 
functional categories were specifically and transiently 
up-regulated in the early germination phase. In contrast, 
the late or post-germination up-regulated genes were 
under-represented in Regulation of Transcription cate- 
gory while late germination down-regulated genes were 
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Figure 7 Over-Representation of Genes Regulated Differentially 

by Each Germination Phase in Master Bins. The representation 

analysis was conducted for the genes differentially regulated by 

each of the three germination phases. Log3 fold change values and 

a Hvu_Affy.m02 mapping file was used in the PageMan analysis. 

Fisher's exact test and an ORA Cutoff value of 1 were used. A false 

color scale of 3 was used to indicate the statistic Z value. Blue 

indicates significance in over-representation while red denotes 

significance in under-representation. Over- and under-represented 

master bins are annotated with functionalities on the left and the 

germination phase and regulation patterns on the top. 
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Figure 8 Over-Representation of Genes Regulated Differentially by Each Germination Phase in All Functional Bins. The representation 
analysis was conducted for the genes differentially regulated by each of the three germination phases. Log3 fold change values and a Hvu_Affy. 
m02 mapping file was used in the PageMan analysis. Fisher's exact test and an ORA Cutoff value of 1 were used. A false color scale of 3 was 
used to indicate the statistic Z value. Blue indicates significance in over-representation while red denotes significance in under-representation. 
Over- and under-represented bins are annotated with functionalities on the right and the germination phase and regulation patterns on the top. 
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over-represented in some regulatory and signaling func- 
tional categories. The specifically and transiently up-reg- 
ulation of the regulatory and signaling functionalities in 
the early germination phase strongly suggested that a 
major molecular event in early germination phase is to 
transcriptionally induce genes encoding regulatory and 
signaling components, and therefore to initiate a variety 
of transcriptional regulatory cascades and signaling 
pathways involved in germination and seedling growth. 

The bin of Cell Wall and its sub-functional bins of Cell 
Wall Precursor Synthesis, Cellulose Synthesis and Cell 
Wall Modification were also up-regulated in the early ger- 
mination phase. A total of 30 cell wall related genes were 
differentially regulated during early germination; and 
twenty-seven of them were up-regulated (See Additional 
file 4). The biased up-regulation of cell wall metabolism is 
consistent with the proposed roles of cell wall synthesis, 
degradation and modification in seed germination. Endo- 
sperm weakening is considered a major feature associated 
with endosperm rupture by the expanding radicle, a major 
regulation point for germination potential [59,63]. In sev- 
eral species endosperm weakening has been associated 
with the induction of cell wall remodeling enzymes. They 
include endo-beta mannanase, beta-l,3-glucanases, expan- 
sins, xyloglucan endotransglycosylase, pectin methylester- 
ase, ploygalacturonase and arabinogalactans [59]. 
Transcripts encoding all of the enzymes except endo-beta 
mannanase were observed in the set of the early germina- 
tion up-regulated transcripts. Barrero, et.al showed that a 
number of cell wall degradation related genes are prefer- 
entially expressed in after-ripening barley coleorhiza for 
both 8 and 18 hours of hydration, and are proposed to be 
associated with breaking seed dormancy [34]. Interestingly, 
the Arabidopsis ortholog of endo-beta mannanase was also 
down-regulated at three hours after imbibition. Fourteen 
of the barley cell wall related genes have Arabidopsis 
orthologous genes that are differentially regulated by more 
than two folds. All of those fourteen pairs of orthologs 
were up-regulated in response to early germination in 
both barley and Arabidopsis, strongly supporting that 
those cell wall genes have functional significance in seed 
germination. It was shown that a GATA zinc finger tran- 
scription factor functions as a positive regulator of germi- 
nation, and is required to facilitate endosperm rupture in 
Arabidopsis [64], Interestingly, three probe-sets (Conti- 
g3743_at, 17684_at and 4186_at) were annotated as 
GATA zinc finger transcription factors and up-regulated 
by more than five folds within the early phase of barley 
germination. Thus, it is likely that early germination pro- 
cess turns on the transcriptional regulatory pathway 
underlying cell wall metabolism activity to weaken coleor- 
hiza and facilitate root emergence. 

The early germination down-regulated genes were 
over-represented in the Stress bin. In this bin, only the 



sub-functional bins of the Heat and Light stresses in 
Abiotic Stress were preferentially down- regulated. Inter- 
estingly, the sub- bins of Light Signaling and Heat-shock 
Transcription Factors were also the only two sub-func- 
tional bins of Signaling and Regulation of Transcription 
over-represented in early germination down-regulated 
genes. The coordinated down-regulation of both heat 
and light-stress pathways and their corresponding tran- 
scription factor and signaling genes in early germination 
phase suggests that biological networks underlying heat- 
and light-stress response were suppressed in early ger- 
mination phase and that down-regulation of those tran- 
scription factors and signaling component genes are 
likely to lead to suppression of the biological networks. 

No significant changes in the morphology, desiccation 
resistance, or amylase activity other than water uptake 
were observed in germinating grains within the first 
nine hours of the early germination phase. However, 
transcript accumulation of 730 genes changed signifi- 
cantly during early germination phase; and germinating 
grains already activated their transcriptional machinery 
and reprogrammed their transcriptional expression to 
synthesize and degrade a specific set of transcripts. 
Interestingly, the studies did not observe that genes 
encoding mRNA synthesis and degradation machinery 
proteins were preferentially up-regulated in the early 
germination phase. It raises the possibility that the ger- 
minating grains utilize RNA synthesis and degradation 
machinery preserved in the mature seeds to support the 
RNA metabolism in the early germination phase. 
Up-regulation of Metabolic Pathways and Chromatin 
Structure in Late-germination Phase 

Over the following nine-hour late germination phase from 
S3 to S4 stages, dramatic morphological and physiological 
changes occurred. Coleorhiza emerged from germinating 
grains; and desiccation resistance decreased significantly. 
However, water-uptake rate reached its lowest level. A 
total of 1295 genes were differentially regulated over the 
late germination process. The bins of Amino Acid Meta- 
bolism, Nucleotide Metabolism, Cell and DNA were pre- 
ferentially up-regulated in the process. In Amino Acid 
Metabolism, the Amino Acid Synthesis pathways, includ- 
ing Methionine and Aromatic Amino Acid Synthesis path- 
ways, but not the Amino Acid Degradation pathways, were 
preferentially up-regulated. Serine Protease in the protein 
degradation pathway, Protein Synthesis Initiation, Major 
Intrinsic Protein Transport and Secretory Pathways in 
Protein Targeting, and Cell Wall Degradation were prefer- 
entially up-regulated in late germination process. In addi- 
tion, the bin of Cell, which includes genes related to cell 
division, cell cycle and vesicle transport, were preferen- 
tially up- regulated in the late germination phase. 

It is striking that the sub-functional bins of Chromatin 
Structure and Histone family were over-represented in 
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late germination up-regulated genes with a very high Z 
value of greater than 18 (See Table 2). A total of 126 
genes encoded histone proteins H2A, H2B, H3 and H4. 
All of the histone genes were up-regulated in the late 
germination phase. In contrast, the functional categories 
of Chromatin Structure and Histone were under-repre- 
sented in both up-regulated and down-regulated genes 
in early and post-germination phases as well as in 
down-regulated genes in late germination phase. Thus, 
histone genes and chromatin structure related genes 
were specifically and preferentially up-regulated in the 
late germination phase. Histone modification and chro- 
matin remodeling play important roles in reprogram- 
ming transcriptional programs, and have been evidenced 
to play regulatory roles in the seed dormancy and ger- 
mination [65-67]. It was shown that the mutation of his- 
tone monoubiquitination genes in Arabidopsis reduced 
ubiquitinated forms of histone H2B, and altered expres- 
sion levels for several dormancy-related genes [66]. A 
transient histone deacetylation event occurs during seed 
germination one day after imbibition, and serves as a 
key developmental signal that affects the repression of a 
number of histone deacetylase regulated genes [67]. 
Gene expression patterns changed very quickly if there 
is no functional constrains after gene duplication [68]. 
The extremely biased up-regulation of so many histone 
genes provide strong evolutionary evidence that they 
might play roles in chromatin remodeling and repro- 
gramming the transcriptional regulatory program for 
seeds to switch from germination to post-germination 
seedling growth. Thus, over-representation analysis data 
suggests that germinating grains activate many meta- 
bolic pathways in late germination phase to produce 
amino acid and nucleotides, establish and maintain pro- 
tein synthesis and transporting machinery, and remodel 
chromatin structure to support the cell division and 
expansion occurring in late and early post-germinations. 

The late germination down-regulated genes were over- 
represented in the functional bins of Fermentation, 
Stress, RNA and Development (Figure 7). Some of the 
sub-functional bins, such as abiotic stress and LEA, 
were down- regulated in both late and post-germination 
phases. Many of the down-regulated functionalities such 
as LEA, Abiotic Stress and ABA Metabolism are highly 



expressed during seed maturation and involved in seed 
maturation and stress tolerance such as desiccation tol- 
erance in mature grains [33,59,69]. Down-regulation of 
those functionalities also co-occurred with loss of desic- 
cation resistance in the late germination period, and is 
likely to contribute to loss of desiccation tolerance dur- 
ing seed germination. 

Up-regulation of Photosynthesis, Degradation, Secondary 
Metabolic Pathways in Post-germination Phase 

Post-germination phase represented the fifty-three hours 
of seedling growth from the S3 to S5 stages. The bins of 
Photosynthesis, Major CHO Metabolism, Secondary 
Metabolism, Tetrapyrrole Synthesis were preferentially 
up-regulated while functional categories of Stress and 
Development were down-regulated (Figure 7). 

In the Photosynthesis pathway, post-germination up- 
regulated genes are over-represented in the sub-bins of 
the Light reaction pathways, the Photosystems I and II, 
Photorespiration and Rubisco Small Subunit Family 
(Figure 8). Figure 9 shows the differentially regulated 
genes in Light Reactions, Calvin Cycle and Photorespira- 
tion pathways in all three germination phases. A total of 
thirty-five genes were differentially regulated in the 
post-germination phase. Thirty-three of the 35 genes 
were up-regulated in the post-germination phase. In 
contrast, early and late germination up-regulated genes 
and early germination down-regulated genes were 
under-represented in the Photosynthesis pathways. 
Thus, the expression of the genes in the photosynthesis 
pathways changed little in the other germination phases, 
and the photosynthesis pathways were specifically up- 
regulated in the post-germination phase (Figure 8). 
Interestingly, genes encoding chloroplast/mito ribosomal 
proteins in protein synthesis pathways were also over- 
represented in post-germination up-regulated genes and 
under-represented in early and late germination up- 
regulated genes (Figure 8). Thus, the photosynthesis 
pathways and chloroplast/mito protein synthesis 
machinery in chloroplast were specifically and coordi- 
nate^ up-regulated in the post-germination phase and 
are likely to support the transition of seed germination 
from heterotrophic growth to photo-autotrophic growth. 
All the samples examined in the studies were germi- 
nated in dark. The specific up-regulation of 
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Figure 9 Differentially Regulated Genes in the Photosynthesis Pathways over Each of Three Germination Phases. The Light Reactions, 
Calvin Cycle and Photorespiration pathways and the differentially regulated genes in those pathways are displayed using the MapMan tool. Each 
differentially regulated gene is displayed as a colored square. Fold changes of their transcript accumulation in each germination process were 
log3 transformed and displayed as a color scale from red (down-regulated) to blue (up-regulated) in the squares. The members of the same 
gene family are grouped and displayed together. Differentially regulated genes in early germination (left, top), late germination (right, top) and 
post-germination are displayed separately. 

V J 



photosynthesis pathways and chloroplast/mito ribosomal 
proteins over post-germination phase suggests that 
plants already develop a machinery or potential in the 
post-germination phase for light response and photo- 
synthesis even in the dark to support their autotrophic 
growth. 

In addition, post-germination up-regulated genes were 
over-represented in many degradation pathways and 
their related transfer protein families. Those include 
Starch Degradation and Starch Cleavage in Major CHO 
metabolism, lipid transfer proteins, beta-oxidation in 
lipid degradation, a variety of amino acid degradation 
pathways, nucleotide degradation, gluco-, galacto -and 
mannosidases, beta 1,3 glucan hydrolyases, nitrilase, acid 
and other phosphatases, GDSL_motif lipase and ribonu- 
cleases. It is consistent to patterns of amylase activity 
over course of seed germination. Up-regulation of indivi- 
dual or limited number of hydrolytic enzyme activities 
and genes during seed germination were documented in 
many publications [24]. However, the over-representation 



analysis of post-germination regulated genes provides 
strong evidence that those seed storage mobilization 
pathways were preferentially up-regulated during post- 
germination at a system level. Interestingly, no protein 
degradation pathway was up-regulated over the post-ger- 
mination phase. Serine Proteases in protein degradation 
pathways were preferentially down-regulated while 
sucrose syntheses, glutamine synthase in N-Metabolism 
and central amino acid synthesis were up-regulated in 
post-germination phase. Thus, those post-germination 
up-regulated pathways and gene families may provide 
additional amino acid and sucrose resource for seedling 
growth. 

Post-germination down-regulated genes were over- 
represented in the bins and sub-bins of Storage proteins, 
LEA and inhibitor proteins, unspecified seed proteins, 
serine protease, raffinose synthesis family, HRGP cell 
wall proteins, TAG synthesis in lipid metabolism, protein 
synthesis, mitochondria targeting pathways, and stress. 
Some of these, such as LEA and stress related 
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functionalities were also preferentially down-regulated in 
late germination. Some of the functional categories, such 
as LEA, TAG synthesis and stresses-related bins were 
highly expressed during seed development and matura- 
tion. It is likely that the mRNA species in the pathways 
are degraded due to decreased needs for their functional- 
ities in seed germination, which is consistent with the 
observation that some seed storage protein transcripts 
are degraded over the course of germination [31]. 

The representation analysis clearly shows that the 
transcriptional program underlying post-germination 
activates various mobilization pathways to degrade the 
storage reserve accumulated in mature grains and to 
meet increasing demands of vigorous seedling growth 
for energy and nutrients. Meanwhile, the growing seed- 
lings also begin to turn on the transcriptional expression 
of photosynthesis pathways to switch from heterotrophic 
growth to photoautotrophic growth. Furthermore, tran- 
scriptional expression of soybean storage proteins, TAG 
synthesis and stress tolerance, which are highly 
expressed over the seed development and maturation, 
were suppressed over the post-germination phase to 
conserve the energy and nutrient for seedling growth. 
Differential Expression Patterns of Biological Pathways 
Interestingly, the up-regulated functionalities were over- 
represented in only one of three germination phases. 
None of them were preferentially up-regulated in more 
than two germination phases, suggesting that each ger- 
mination phase transcriptionally induces a distinct set of 
biological functionalities over the course of germination. 
However, several functional categories related to stres- 
ses, development and LEA were down regulated over 
more than two phases. Most of the mRNA species in 
those functional categories were synthesized over the 
seed maturation and often highly accumulated in mature 
grains. They are likely to be degraded progressively in a 
less strictly regulated manner over the course of seed 
germination. However, no bin or sub-bin was over- 
represented in both up- regulated genes and down-regu- 
lated genes in the same germination phase. However, a 
number of functional categories were over-represented 
in up (down) -regulated genes, but under-represented in 
down (up)-regulated genes in the same germination 
phase. This suggests that plants have evolved a sophisti- 
cated and well regulated transcriptional mechanism to 
suppress antagonistic regulation of the same functional 
categories in the same germination phase to efficiently 
use the energy to support seed germination and seedling 
growth. 

Transcriptional Regulatory Programs Underlying GA and 
ABA Regulation of Seed Germination 

We previously examined transcriptomes of isolated ger- 
minating barley aleurone treated with GA and ABA 



respectively, and identified 1328 GA-responsive genes 
and 206 ABA-responsive genes in their mRNA accumu- 
lation [20]. Comparing the GA or ABA responsive genes 
with 4493 germination responsive genes showed that 
46% of the GA responsive genes and 57% of the ABA 
responsive genes were also differentially regulated dur- 
ing germination (Figure 10A and 10B), suggesting that a 
large portion of GA or ABA responsive genes identified 
in isolated aleurone cells are likely involved in seed 
germination. 

A total of 610 genes were differentially regulated by 
both the GA treatment and germination. Of the 610 
genes, 490 (80%) genes showed a coordinate response to 
GA and germination (See Additional file 5), while only 
128 genes showed an antagonistic response to GA treat- 
ment and germination (Figure 10A). The number of 
coordinately regulated genes was 3.8 times that of antag- 
onistic regulated genes. The preferentially coordinate 
response of the genes to GA and germination provides 
strong evidence at a systems level for the hypothesis 
that GA enhances seed germination and seedling growth 
[4,5,12]. The coordinately regulated genes are potentially 
involved in GA regulation of germination and seedling 
growth processes. 

Representation analysis revealed that the major CHO, 
cell wall and protein degradation pathways, which also 
included starch degradation, starch cleavage, cell wall 
degradation, cysteine protease and ubiquitin degradation 
pathways, were preferentially up-regulated by both GA 
and germination. However, TAG synthesis, Aspartate 
family degradation, metal handling, Hormone metabo- 
lism, ABA induced genes, Short chain dehydrogenase/ 
reductase (SDR), and Late embryogenesis abundant 
genes were preferentially down-regulated by both GA 
and germination. Thus, it is likely that GA enhances 
seed germination and seedling growth partly through 
inducing mobilization of cell wall, starch and protein, 
and suppressing production of many seed maturation 
proteins such as LEA, TAG synthesis and ABA activated 
functions. Interestingly, the studies did not observed 
that the genes antagonistically regulated by GA and ger- 
mination were over-represented significantly in any 
functional bins or sub-bins (Figure 11). 

Out of 206 ABA responsive genes, 118 were also dif- 
ferentially regulated by germination. 95 of the 118 genes 
(80%) showed an antagonistic response to ABA treat- 
ment and germination (Figure 10B). Sixty-four genes 
were up-regulated by ABA but down-regulated by ger- 
mination. Many of them encode embryogenesis abun- 
dant proteins and stress-related proteins. For example, 
11 late embryogenesis abundant protein genes, 3 dehy- 
drin genes (Contigl721_at, Contigl718_s_at, Conti- 
gl709_at) and three Glutathione-S-transferase genes 
(Contigl4304_at, Contig2248_at, Contig2975_s_at) (See 
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10A 



490 Up-Up or Down-Down 



1328 GA 




4492 Germination 



128 Up-Down or Down-Up 



95 Up-Down or Down-Up 



10B 



206 ABA 




4492 Germination 



24 Up-Up or Down-Down 



4492 Germination 



IOC 



1328 GA 




206 ABA 



43 GA/ABA antagonistic, ABA/germination antagonistic and 
GA/germination coordinated probe-sets; 13 probe-sets for the rest. 

Figure 10 Venn Diagram of the Genes Differentially Regulated by Germination, GA or ABA. Germination, GA and ABA differentially 
regulated genes were compared and displayed as Venn diagrams. 10A compares GA and germination; 10B compares ABA and germination; and 
10C shows the comparison of GA, ABA and germination. The number of probe-sets and their expression patterns are shown. 
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Figure 11 Over-Representation of the Genes Differentially Regulated by Both GA and Germination. The representation analysis was 
conducted on the 490 genes coordinately regulated by both GA and germination. Fold changes in their transcript accumulation in response to 
GA treatment was used in this analysis. Fisher's exact test and ORA Cutoff value of 1 were used. A false color scale of 3 indicates the statistic Z 
value. Blue indicates significance in over-representation while red denotes the significance in under-representation. Over- and under-represented 
bins are displayed using PageMan Tool and annotated with functional categories on the right and the expression patterns on the top. 



Additional file 6). The genes down-regulated by ABA 
and up-regulated by germination included alpha amy- 
lase, 3 beta-glucanase, beta-xylanase and cysteine protei- 
nases. It has been well established that ABA promote 
establishing and maintaining seed dormancy, and inhi- 
biting seed germination [13,18]. The studies provide sys- 
tems evidences supporting that ABA suppresses 
germination partly through inhibiting storage reserve 
mobilization and enhancing expression of maturation 
genes. 

A total of 56 germination responsive genes were dif- 
ferentially regulated by GA and ABA (See Figure 10C 
and See Additional file 7). Forty- three out of the 56 ger- 
mination responsive genes were regulated coordinately 
by GA, and antagonistically by ABA. Seventeen of the 
genes were up-regulated by both GA and germination, 
but down-regulated by ABA. Twenty six genes were 



down-regulated by both GA and germination, but up- 
regulated by ABA. Many of the genes encoded embryo- 
genesis abundant and stress-related proteins, which 
include 10 LEA proteins, one GST and one dehydrin. 
Four regulatory genes encoding two protein phosphatase 
2C, a WRKY transcription factor and an AP2-domain 
DNA-binding protein were identified in the set of genes. 
It has been proposed that GA and ABA play antagonis- 
tic roles in regulating seed maturation and germination 
[4,5]. Those genes are likely to be important part of the 
pathway mediating the antagonistic interaction of GA 
and ABA signaling pathways in regulating seed germina- 
tion. The four transcription factor and signaling genes 
may play a regulatory role in the networks. It should be 
noted that those GA and ABA responsive genes are 
identified in barley aleurone. However, GA and ABA 
regulation of seed germination and seedling growth 
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could occur at multiple tissues including embryo tissues. 
It remains to be determined if those genes are differen- 
tially regulated by GA and ABA in those tissues as in 
barley aleurone. 

Conclusions 

The studies carefully examined water content, loss of 
desiccation tolerance, amylase activity and morphology 
of germinating barley, and selected six distinct develop- 
mental stages for transcriptome analysis based on the 
multi-germination characteristics and relative timing of 
germination. The studies developed a model depicting 
transcriptional regulatory program underlying barley 
germination at gene, pathway and systems levels. Prior 
to grain imbibition, mature barley grains already accu- 
mulate a large number of transcripts, which are synthe- 
sized during seed development and maturation and 
preserved in mature grains. Although accumulation of 
some of those transcripts decreases over the course of 
barley germination, a significant number of the tran- 
scripts are likely to remain in germinating barley and 
participate in seed germination and seedling growth. 
Comparing highly abundant transcripts in barley and 
Arabidopsis dry seeds showed that those barley and Ara- 
bidopsis transcripts in dry seeds are highly conserved, 
and suggested the ancient origins of those highly abun- 
dant seed transcripts and their functional significance in 
germination. Upon grain imbibition, a new transcrip- 
tional regulatory program is activated quickly, which 
could occur as early as within the first three hours of 
imbibition. Hierarchal clustering of the transcriptomes 
of germinating barley at each developmental stage 
reveals that the new transcriptional regulatory program 
is composed of three distinct phases, early, late, and 
post-germination phases. Early germination phase repre- 
sents the first nine-hour germination and preferentially 
induces genes encoding regulatory components includ- 
ing transcription factors, signaling components and 
post-translational modification proteins. Those regula- 
tory genes are likely to activate a variety of transcrip- 
tional regulatory cascades and signaling transduction 
pathways in seed germination and seedling growth. In 
addition, cell wall synthesis and modification pathway 
genes are also preferentially up-regulated within the 
early germination phase, which may function to loosen 
cell walls for subsequent cell expansion and division, 
and radicle protrusion. Within the following 9 hours of 
late germination phase, genes encoding many metabolic 
pathway enzymes and cellular components including 
amino acid and nucleotide synthesis, protein degrada- 
tion, chromatin remodeling and cell division pathways 
are preferentially up-regulated to provide nutrient and 
cellular components for cell division and elongation. In 
addition, a transcriptional switch occurs in the late 



germination phase, and correlates with the developmen- 
tal transition from germination to seedling growth. 
Post-germination phase mainly represents seedling 
growth process after coleorhiza emergence from the ger- 
minating grains. As expected, seed reserve mobilization 
and photosynthesis pathway genes are preferentially 
induced in the phase to mobilize seed storage reserve in 
mature grains, and acquire autotrophic growth ability to 
meet the increasing demand for energy and nutrients in 
seedling growth. In addition, transcriptional expression 
of secondary metabolism and tetrapyrrole synthesis 
pathway genes are up-regulated during post-germination 
phase. However, many genes encoding stress related 
protein, LEA genes and seed storage proteins, which are 
highly expressed during seed development and matura- 
tion, are transcriptionally suppressed over the course of 
barley germination to conserve the energy and nutrients 
for seed germination and seedling growth. Although 
some of those metabolic pathways have been previously 
proposed or logically assumed to play roles in germina- 
tion, the studies illustrated that transcriptional expres- 
sion of those pathways are differentially regulated over 
the course of germination at a systems level, and pro- 
vided additional evidences for their roles in germination. 
The studies also newly discovered a number of pathways 
and genes that are differentially regulated over the 
course of germination at a transcriptional level, and sug- 
gested their functional involvement in germination. A 
great number of hypotheses also have been developed in 
the studies for future validation. In addition, the studies 
identified a set of genes encoding regulatory compo- 
nents that were transiently up-regulated as early as 3 
hours of imbibition. Their transient and up-regulated 
expression patterns at such an early stage of germina- 
tion suggests that they may play key regulatory func- 
tions in seed germination, and worth further 
investigating their functions in seed germination. 

The studies also compared GA and ABA responsive 
genes with genes differentially regulated by barley ger- 
mination, and identified three sets of germination 
responsive genes that coordinately respond to GA, 
antagonistically respond to ABA, and coordinately 
respond to GA but antagonistically respond to ABA. 
Those genes are likely to be important components in 
the transcriptional regulatory networks that GA 
enhances germination, ABA suppresses germination, 
and GA and ABA interact antagonistically in regulating 
germination. 

Overall, the studies establish a standard transcriptome 
reference platform for barley germination and enable 
seed biologists to integrate the transcriptome data with 
a variety of -omics and other biological data to illustrate 
biological networks underlying barley germination. The 
studies also developed a model depicting the 
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transcriptional regulatory programs underlying seed ger- 
mination and germination related biological pathways, 
and GA and ABA regulation of seed germination at 
gene, pathway and systems levels. 

Methods 

Plant Growth and Harvest 

Plump and healthy barley grains (Hordeum vulgare L. 
cv. Morex L.) were imbibed in water for three hours at 
22°C with three changes of water and were then germi- 
nated on water-saturated germination pack in the dark 
at 22°C. Twenty grains were planted in each 15 cm dia- 
meter Petri-dish and spaced evenly to reduce the varia- 
tion of seed germination caused by grain density. Over 
98% of grains germinated at these conditions, but indivi- 
dual grains did not germinate uniformly. The germinat- 
ing grains/developing seedlings with typical 
morphologies at 0 (dry), 3, 9, 18, 33, and 71 hours of 
germination were harvested and pooled for determina- 
tion of water content, alpha amylase activity, loss of 
desiccation resistance and RNA purification. The typical 
morphology of the seedlings at 18, 33, and 71 hours of 
germination was identical to that of S3, S4 and S5 
stages, (Figure 1A). The grains/seedlings harvested for 
RNA extraction and determination of alpha amylase 
activity were pooled together and immediately frozen in 
liquid nitrogen and stored at -80°C. Each replication of 
water content, amylase activity, desiccation resistance 
and microarray assays represented an independent ger- 
mination experiment. Three independent germination 
experiments were conducted for each time point. Each 
replication represents a pool of germinating barleys that 
were carefully selected from an independent germina- 
tion experiment based on the relative time point and 
morphology defined for the given stage to reduce the 
heterogeneity of the pooled barley tissues potentially 
caused by their different germination rate. 

Alpha-Amylase Assay 

Alpha amylase activity was measured using a DNSA 
assay [70]. One-half gram of tissue was homogenized 
using a Bead-Beater at maximum speed for 2 min in 
500 ml of phosphate buffer (20 mM Na 2 HP0 4 , 10 mM 
NaCl, pH 6.9). The sample was then centrifuged for 10 
min at 13,000 rpm. The aqueous phase was incubated at 
69°C for 15 min to inactivate beta-amylase. After centri- 
fugation for 10 min at 13,000 rpm, 10 mL of solution 
was added to a phosphate buffer with 0.5% starch and 
incubated at 30°C for 30 min. Then, an equal amount of 
DNSA reagent (0.25 mM NaO, 1% 3,5-dinitrosalicylic 
acid, 30% NaK tartrate) was added to the reaction and 
incubated for 15 min at 100°C. After cooling the sample 
to 22°C and centrifuging for 5 min at 13,000 rpm, 200 
ml of supernatant was added to microplates to measure 



the OD at the wavelength of 547 nm. Maltose was used 
as the standard to calculate enzyme activity. Dry grains 
were used as a control for amylase activity. 

Desiccation Resistance Assay 

For the desiccation resistance assay, the grains/seedlings 
with the representative morphology at each given stage 
(Figure 1A) were harvested and dried aerially at room 
temperature for two weeks. The dehydrated grains/seed- 
lings were re-germinated at the same condition as 
described in the germination experiment. The desicca- 
tion resistance of the grains at each stage was deter- 
mined as percentage of survival, defined as the 
percentage of the dehydrated grains/seedlings that can 
revive within three days of germination. Sample size for 
each desiccation resistance assay ranged from 23 grains 
to 56 grains. Three independent replications were con- 
ducted for each time point. 

RNA Purification 

Total RNA extraction was conducted as described by 
Chen and An with minor modifications [20]. Two grams 
of plant tissue was ground in liquid nitrogen followed by 
the addition of 10 ml extract buffer (4% ^-aminosalicylic 
disodium, 1% 1,5-napthalenedisulfonic acid) and 10 ml 
phenol. The mixture was inverted several times, and 
polytroned for 45 seconds after the addition of 10 ml 
chloroform. After centrifuging, the aqueous phase was 
transferred into a new tube, to which 60 ml of 10% Cal- 
coflur White was added [70]. The mixture was mixed 
thoroughly and centrifuged for 15 min at 4°C, 12,000 
rpm. Precipitate RNA from the supernatant was formed 
by using 1/10 volume of 3M NaO Ac, and 2 volume of 
100% ethanol. After centrifuging, the RNA pellet was 
first dissolved in 8 ml water, then 5 ml of 8M LiCl was 
added, and the tubes were left on ice overnight. After 
centrifuging, the resulting RNA pellet was dissolved in 
water. The RNA quality and quantity was checked by 
using Nano-Drop (Nano-Drop, Wilmington, DE) and 
Aglient 2100 Bioanalyzer (Aglient, Palo Alto, CA). 

GeneChip Array Assay 

Preparations of cDNA and biotin-labeled cRNA were 
performed and analyzed as recommended by Affymetrix, 
Inc. (Santa Clara, CA). According to the manufacturers 
protocol, 7.5 to 15 mg of total RNA was used in a 
reverse transcription reaction to generate first-strand 
cDNA using Reverse transcriptase Superscript II (Invi- 
trogen, Carlsbad, CA). After second-strand synthesis, 
double-strand cDNAs were used in an in vitro transcrip- 
tion reaction to generate biotinylated cRNA. Ten mg of 
fragmented cRNA was used for each hybridization. 
Staining and scanning of the hybridized GeneChips 
were performed as described by Chen and An [20] and 
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the manufacturer's recommended protocols (Affymetrix, 
Inc., Santa Clara, CA). The .CEL files from all samples 
are available in GEO (http://www.ncbi.nlm.nih.gov/geo/ 
info/linking. html) with series record number of 
GSE23595 and access numbers of GSM578686 to 
GSM578703. 

Data acquisition and analysis 

Affymetrix GeneChip Microarray Suite version 5.0 soft- 
ware (MAS 5.0) was used to assign the presence and 
absence calls of each probe set for each GeneChip with a 
P value of 0.05. The data files containing the probe set 
intensities (.eel files) were used for background correction 
and normalization by an improved log2 scale RMA proce- 
dure, GC-RMA, provided in GeneSpring Suite 7.2. Within 
each array, a further "per gene normalize the median" 
(with cutoff 0.01) was applied to the pre-normalized data 
using GC-RMA provided in the GeneSpring 7.2 software. 
The probe sets with absence calls across all chips were 
removed from further analysis. A probe set with present 
(or absent) calls in two of the three replicates were 
assigned as present (or absence) call for the treatment. We 
used the GC-RMA approach to convert probe level data 
to expression measurement in the microarray experi- 
ments. Compared with the algorithms used in Microarray 
Suites (MAS 5.0), this approach adjusts background on 
the raw intensity value scale, then uses quartile normaliza- 
tion to remove systemic variations, and summarizes log2 
of the normalized background adjusted PM values to esti- 
mate expression level measurements based on a linear 
additive model [39] . One-way ANO VA was used to iden- 
tify genes that were differentially expressed at any two 
time points during seed germination with a False Discov- 
ery Rate (FDR) of 0.05. The Parametric Test, Variances 
Assumed Equal Option, Benjamini and Hochberg multiple 
testing corrections were used in the one-way ANOVA 
analysis. The harvEST:Barley (version 1.35; http://harvest. 
ucr.edu/), Munich Information Center for Protein 
Sequence [71] and Universal Protein Resource [72], were 
used to conduct gene functional annotation in addition to 
manual editing. An E score of 1 x E" 20 of BLAST-X 
between a barley sequence and Arabidopsis sequences in 
HarvEST: Barley was used as a cutoff.. The Arabidopsis 
genes homologous to a given barley gene with the lowest 
e-value were defined as putative Arabidopsis orthologs. 
MapMan (version 1.4.3) and PageMan (Version 0.12) were 
used for over- representation functional analysis [62,73]. 

Additional material 



Additional file 1: Genes differentially regulated between any stages 
vs. SO and between adjacent stages. 

Additional file 2: Summary of Differentially Regulated Genes over 
the Course of Barley Germination. 



Additional file 3: List of Probe-sets Differentially Regulated Based 
on Germination Phases. 

Additional file 4: Early Germination Regulated Cell Wall Genes. 

Additional file 5: Genes Differentially Regulated by Both GA and 
Germination Genes. 

Additional file 6: Genes Differentially Regulated by Both ABA and 
Germination. 

Additional file 7: Genes Differentially Regulated by GA, ABA and 
Germination. 
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